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We report the low-frequency and tunneling studies of yttrium hexaboride single crystal. Ac 
susceptibility at frequencies 10 - 1500 Hz has been measured in parallel to the crystal surface DC 
fields, Ho. We found that in the DC field Ho > H C 2 DC magnetic moment completely disappears 
while the ac response exhibited the presence of superconductivity at the surface. Increasing of 
the DC field from H C 2 revealed the enlarging of losses with a maximum in the field between H C 2 
and H C 3. Losses at the maximum were considerably larger than in the mixed and in the normal 
states. The value of the DC field, where loss peak was observed, depends on the amplitude and 
frequency of the ac field. Close to T c this peak shifts below H C 2 which showed the coexistence 
of surface superconducting states and Abrikosov vortices. We observed a logarithmic frequency 
dependence of the in-phase component of the susceptibility. Such frequency dispersion of the in- 
phase component resembles the response of spin-glass systems, but the out-of-phase component 
also exhibited frequency dispersion that is not a known feature of the classic spin-glass response. 
Analysis of the experimental data with Kramers-Kronig relations showed the possible existence of 
the loss peak at very low frequencies (< 5 Hz). We found that the amplitude of the third harmonic 
was not a cubic function of the ac amplitude even at considerably weak ac fields. This does not 
leave any room for treating the nonlinear effects on the basis of perturbation theory. We show 
that the conception of surface vortices or surface critical currents could not adequately describe 
the existing experimental data. Consideration of a model of slow relaxing nonequilibrium order 
parameter permits one to explain the partial shielding and losses of weak ac field for Ho > H C 2- 

PACS numbers: 74.25.Nf, 74.25.Op, 74.70.Ad 
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I. INTRODUCTION 



There is a growing interest in exploring physical prop- 
erties of materials involving boron-cluster compounds be- 
cause of a wide variety of applications [l[ . Due to sp 2 hy- 
bridization of valence electrons, large coordination num- 
ber and short covalent radius, boron atoms prefer to form 
strong directional bonds with various elements. A large 
number of experimental and theoretical studies are con- 
centrated on the families of compact B12 icosahedrons 
and B6 octahedrons with a large diversity of electrical 
and magnetic characteristics. The highest critical tem- 
peratures of the transition to the superconducting state 
in MBg and MB12 compounds were found in YBq with 
T c < 8.4 K and ZrBi 2 with T c w 6.0 K @|. Both materials 
have a highly symmetrical crystal structure (CaBgtype 
for YB 6 and UB12 type for ZrBi 2 ) that can be described 
as boron cages in which yttrium or zirconium atoms de- 
velop large vibrational amplitudes with an Einstein-like 
(nearly dispersionless) lattice mode. In spite of some 
common features, the two crystals have a few distinct 
physical characteristics: (i) while YBg is a classical type- 
II superconductor [3j], ZrBi2 (at least, for temperatures 
above 4.5 K) may be regarded as a textbook example of 



type-I superconductor [4| ; (ii) while the superconducting 
properties are enhanced at the ZrBi2 surface [4J, they 
are suppressed in a YBg surface (see our tunneling data 
below). Therefore, ZrBi 2 and YBg samples may serve 
as model systems for investigating surface-related super- 
conducting effects. 

Nucleation of a superconducting phase in a thin sur- 
face sheath, when the DC magnetic field, Hq, paral- 
lel to the sample surface decreases, was predicted in 
1963 by Saint- James and de Gennes in their seminal 
work fc 5j. They showed that the nucleation occurs for 
Hq < H C 3 = 2.39kH c , where H c is the thermodynamic 
critical field and k is the Ginzburg-Landau(GL) param- 
eter. Experimental measurements confirm this predic- 
tion @, S, E3, El [3, and it was found that at low 
frequencies a sample in a surface superconducting state 
(SSS) shows ac losses with a peak whose position with 
respect to the DC field depends on the ac amplitude. The 
peak magnitude exceeds the losses observed either in the 
normal state (Hq > H c ^) or in the bulk superconducting 
state(-f/o < H C 2). It was also predicted that the H^/H^ 
ratio, is temperature independent. In contrast, a decrease 
of this ratio was found in the vicinity of T c in several 
experiments [111, [l2| ■ This behavior was associated with 
the distribution of T c at the surface [l3T |. 
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In the last few years the SSS has attracted renewed 
interest from various directions 0, EE EE EE EE El- 
Stochastic resonance phenomena in Nb single-crystal 
were observed in the SSS In [TtJ it was assumed 

that at Ho < H c ^ the sample surface consists of many dis- 
connected superconducting clusters and subsequently the 
percolation transition takes place at H£ 3 = 0.81-ff C 3. The 
paramagnetic Meissner effect is also related to the SSS 
|18j . Voltage noise and surface current fluctuations in Nb 
in the SSS have been investigated [Hj]. SSS' were found 
also in single crystals of ZrB 12 (lq |. In agreement with 
the previous data @ it was demonstrated [l6| that the 
waveform of the surface current in an ac magnetic field 
has a non-sinusoidal character. A simple phenomenolog- 
ical relaxation model provides the good explanation of 
the experimental data for DC fields near H c i only 
The relaxation rate in this model depends on the ac fre- 
quency and decreased with decreasing ui pj|. Detailed 
experimental study of the linear ac response in the SSS 
of single crystals Nb and ZrBi2 was published recently in 
our paper [4| . We showed that ac SSS losses in these ma- 
terials could be considered in the achieved experimental 
accuracy as a linear ones and for several DC fields the 
real part of the ac magnetic susceptibility exhibited a 
logarithmic frequency dependence as for a spin-glass sys- 
tem. 

In spite of the extensive studies, the origin of low fre- 
quency losses in SSS is not clear as yet. The critical state 
model developed for the SSS in [19| implies that if the 
amplitude of the ac field, ho, is smaller than some crit- 
ical value the losses disappear. The authors of Ref. @ 
claimed that the experiment on Pb-2%In alloy confirms 
this prediction. On the other hand, the observed re- 
sponse p3| for an excitation amplitude of 0.01 Oe that is 
considerably smaller than used in Q showed losses in SSS 
in Nb sample at a frequency 10 Hz. Our measurement 
on Nb and ZrB2 single crystals also have shown that the 
out-of-phase part of the ac susceptibility, was finite 
at low excitation level [J] . We consider these results as an 
indication of the inadequacy the critical state model for 
description of the ac response in SSS. If we assume that 
the reason for this discrepancy with experimental data, 
is the small value of the critical surface current, which is 
much smaller than the current amplitudes at the surface, 
then we have to expect a decreasing of the losses approx- 
imately as 1/ho when the amplitude of the applied ac 
field, ho, increases. On the contrary, \\ increases with 
Iiq. The ac investigation of YBg samples has some advan- 
tage due to actually ideal type II magnetization curves 
in this material that permits one to avoid possible dif- 
ficulties in the interpretation of the experimental data. 
The experiment showed that near the transition temper- 
atures SSS exist also in the fields below H C 2. We found 
that some features of nonlinear response took place at 
a very weak ac field with amplitude ~ 0.005 Oe. The 
ac response at the third harmonic of the fundamental 
frequency did not leave any room for the perturbation 
theory. It was proposed that the losses in SSS are due 



to the slow relaxation of the order parameter at the sur- 
face and could not be ascribed to surface vortices. We 
found that for small ho in quasilinear approximation the 
integral equation with power dependent nuclear governed 
the time behavior of the magnetization in ac fields. Some 
features of the ac response resemble the ones of spin-glass 
system but one has to note that SSS present a different 
system with its unique properties. 



II. EXPERIMENTAL DETAILS 

A. Sample preparation 

The yttrium hcxaboride single crystal was grown by 
the inductive floating zone method of a powder sintered 
rod with an optimal composition YBg.ss under 1.2 MPa 
of argon. According to the Y - B phase diagram, com- 
position with the Y:B=1:6 ratio has undergoes peritectic 
melting [2(| and irrespective of the Y/B ratio the YB4 
single crystal with preferential orientation [001] begins 
to grow. After enrichment of the melting zone by boron 
(flux method modification) the yttrium hexaboride single 
crystals with the [100] orientation grow with the compo- 
sition of YB5.79io.02 (ESD). The total impurity concen- 
tration is less than 0.001 % in weight and the obtained 
lattice parameter is 4.1001(4) A in accordance with pub- 
lished data [2l|. These single crystals exhibited a sharp 
superconducting transition with T c w 7.15 K. 



B. DC and ac measurements 

The magnetization curves were measured using a com- 
mercial SQUID magnetometer. In-phase and out-of- 
phase components of the ac susceptibility at the funda- 
mental frequency, and the response at the third harmonic 
were measured using the pick-up coil method [E I22I ]. 
A home-made setup was adapted to the SQUID mag- 
netometer, and the block diagram of the experimental 
setup was published in Ref. [16[. The crystal (10 x 3 x 1 
mm 3 ) was inserted into one of a balanced pair of coils. 
The unbalanced signal and the third harmonic signal as 
a function of the external parameters such as temper- 
ature, DC magnetic field, frequency and amplitude of 
excitation, were measured by a lock-in amplifier. The 
experiment was carried out as follows. The sample was 
cooled down at zero magnetic field (ZFC). Then the DC 
magnetic field, Ho, was applied. The amplitudes and the 
phases at all frequencies of both signals were measured 
in a given Ho (including at zero field). The excitation 
amplitude, ho, was 0.0005 -j- 0.5 Oe. It is assumed that 
in Ho — 0, and at low temperature, the ac susceptibil- 
ity equals to the DC susceptibility in the Meissner state 
with negligible losses. This permits us to find the abso- 
lute values of the in-phase and out-of-phase components 
of the ac susceptibility for all applied DC and ac fields 
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and for all frequencies. Both Hq and ho were parallel to 
the longest sample axis. 



C. Tunneling measurements 

Measurements of the tunneling spectra were carried 
out using a home made scanning tunneling microscope. 
The YBg single crystal was mounted inside the cryogenic 
scanning tunneling microscope and then cooled down to 
4.3 K. The dl/dV vs. V tunneling spectra (proportional 
to the local density of states) were acquired using a con- 
ventional lock-in technique, while momentarily discon- 
necting the feedback loop. 
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III. EXPERIMENTAL RESULTS 

A. Tunneling characteristics 

Direct information about the energy gap value A = 
A(T = 0) at the surface of YB 6 was obtained from the 
tunneling spectra. The ratio 2Aq/T c is a well known 
indicator of the electron-phonon coupling strength [24| . 
Two previous tunneling studies of YBg were performed 
on a single-crystal Q and on thin films [23| and the Aq 
obtained were: 1.22 [f] and 1.24 [H| which yields the 
ratio 2Ao/T c sa 4. In both cases the tunneling con- 
tacts were connected to underlying layers, and hence, 
monitored bulk properties. Therefore those values which 
signified a nearly strong coupling are attributed to the 
bulk characteristics. It was confirmed, in particular by 
Lortz et al. [25j |. who measured the deviation function 
D{T) = H C (T)/H C (G) - (1 - (T/T c ) 2 ) and found that 
the value of 2Ao/T c is slightly above 4.0. Our tun- 
neling spectroscopy results were obtained by STM and 
therefore better reflect the density of states at the sur- 
face. In contrast to our previous measurements on ZrBi2 
single crystals that showed very high spatial homogene- 
ity [H, the superconductivity in the present case ap- 
peared to be degraded on parts of the YBg sample sur- 
face, where nearly featureless tunneling spectra were ob- 
served. In other regions, however, reproducible ratios 
of differential conductances in superconducting and nor- 
mal states (dI/dV) s /(dI/dV) n showing very clearly that 
BCS-like gap structures were acquired, such as presented 
in Fig. [T] (solid line). The spectra were compared with a 
temperature-smeared version of the Dynes formula [26| 
wich takes into account the effect of incoherent scatter- 
ing events by introducing a damping parameter T into the 
conventional BCS expression for a quasiparticle density 
of states 

N S (E) = N N (0)Re[(E-iT)/y/(E - iT f - A*(T)]. (1) 

A very good fit to the experimental data (except for a 
small asymmetry in the normal resistance between neg- 
ative and positive bias, the origin of which is not yet 



FIG. 1: Representative tunneling spectrum of YBg at T 4.3 
K (solid line) together with its fit to the Dynes function (see 
text) shown by a dashed line, with fitting parameters A(T) = 
1.0 meV and V = 0.10 meV. The spectra were normalized to 
the normal tunneling conductance at 5 meV (well above the 
superconducting gap) . 



clear to us) was achieved with A(4.3if) = 1.0 meV and 
r = 0.10 meV. Recalling that the experimental spec- 
trum was acquired at T — 4.3 K, which is about 0.6T C , 
with the BCS A(T) dependence [231 we obtain the zero- 
temperature value A(T = 0) = 1.1 meV. With that, we 
find that 2Ao/T c 3.59, very close to the BCS weak 
coupling value of 3.53. In contrast to ZrBi2, we assume 
that in YBg the electron-phonon strength is suppressed 
at the surface to a weak coupling state. 



B. DC and ac magnetic characteristics 

Fig. [^demonstrates the temperature dependence of the 
sample magnetic moment. In this curve one can see that 
T c w 7.15 K. 

From the hysteresis curve measured at 4.5 K, shown 
in Fig. [31 we are able to evaluate H c = 295 Oe, H C 2 = 
1500 Oe and GL parameter m = H c2 /V2H C = 3.58. Us- 
ing the relation ^jj-\H a =H a2 = l/47rA4(2«;| — 1) on e can 
obtain that k 2 = 3.3, where 0a = 1.16. The temperature 
dependence of H C 2 is shown in the inset of Fig. [3] The 
London penetration depth at T=0, Al(0), can be esti- 
mated by using H c2 (T) near T c , dH c2 /dT w -560 Oe/K, 

VAl(0) = ^g\W\ 0, and A L (0) « 1.4 x lO" 5 cm. 

Fourier analysis of the magnetization, M (t) , un- 
der applied ac and DC fields, H(t) = Hq + 
hocos(ut), yields an expression: M(t) — Mo(H ,h ) + 
Tin ^Xn(Ho, ho)ho exp(-inujt) . In this paper we discuss 
the results for \i an d X3 susceptibilities. The field de- 
pendence of xi(Ho) at T = 4.5 K and ho = 0.05 Oe for 
some frequencies is shown at Fig. 01 One can readily see 
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FIG. 2: (Color online) Temperature dependence of magnetic 
moment after ZFC. 



FIG. 4: (Color online) Magnetic field dependencies of Xi and 
Xi at T — 4.5 K at different frequencies u). 
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FIG. 3: (Color online) Magnetization curve at T — 4.5 K 
after ZFC. Inset: temperature dependence of H C 2- 



that the curves shift toward higher DC fields with fre- 
quency. Decreasing the ac amplitude produces a similar 
effect. The curves shift to the higher field when ho — > 
( Fig. [5]). Similar effects were reported for a Pb-2%In 
sample in Ref. @ . 

The typical magnetic field dependence of the nonlinear 
response, X3, is presented at Fig.[|5]for h — 0.05 Oe and 
various frequencies. When the frequency increases the 
maximum in X3 moves toward larger DC fields as was 
observed for x'i ■ The frequency dispersion is illustrated 
on the Cole-Cole plot, Fig. [7] One can see x'i (panel a) 
and X3 (panel b) as a function of x'i when the frequency 
increases from 15 to 1465 Hz while the DC field was kept 
constant. Each disconnected curve of this figure corre- 
sponds to different DC fields, the values of which are 
indicated in panel (b) . The arrow in panel (b) shows the 
direction of increasing frequency along the curves and 
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FIG. 5: (Color online) Magnetic field dependencies of Xi and 
Xi at T — 4.5 K at different amplitudes of excitation, ho. 



shielding, as well as \xx\- Below H C 2 iirx'i = — 1 (see Fig. 
4). For H > H C 2 both x'i and X3 decrease as the fre- 
quency increases while for Ho close to H c s they increase. 

Fig.[8]shows the field dependence of X3 at u/2tt = 1465 
Hz and various amplitudes of excitation, ho- The third 
harmonic cannot be adequately described in the frame of 
the perturbation theory which predicts that X3 oc h\. For 
example, at Hq/H C 2 — 1.3, X3 depends on ho strongly, 
while at H /H c2 ~ 1.45, X3 1S almost constant (see 
Fig. [5]). We can discuss only the dependence of X3m (de- 
fined as the maximum value of the X3(^o) curve for any 
given frequency) on the ac amplitude h - Fig. [3] demon- 
strates that X3m ~ ft-o 2 in contrast to what the pertur- 
bation theory predictions. 

Below we consider the experimental results obtained 
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FIG. 6: (Color online) Third order susceptibility, X3, versus FIG. 8: (Color online) Third order susceptibility, X3, versus 
reduced magnetic field, H /H c2 , at different frequencies. reduced magnetic field, H /H c2 , at different amplitude of ex- 

citation. 





FIG. 7: (Color online) Panel (a): Cole-Cole plot of the first 
harmonic ac susceptibility. Panel (b): X3 versus x'i- F re ~ 
quency ui and DC field Ho are parameters for these parametric 
curves. The symbols on the both panels are the same. 



FIG. 9: (Color online) Amplitude dependence of the third 
order susceptibility at maximum, X3m, at different frequencies 
(see text). 



at higher temperatures. Fig. [TU] demonstrates the field 
dependence of \i & t frequency u/2ir = 1065 Hz and 
ho = 0.05 Oe at various temperatures. The peak in Xi 
shifts toward H C 2 with temperature and at 7 K this peak 
is located already below H C 2 . One can see in the Fig. [10] 
that for T < 7 K full shielding (xi = — l/47r) is observed 
at low Ho, whereas at 7 K only partial shielding is ob- 
served at low DC field. Also Fig. [TT] shows that in the 
vicinity of T c Xzm lies below H C 2- Because we did not ob- 
served any absorption peak and harmonic signal in the 
mixed state we consider that SSS are responsible for the 
experimental observations at T = 7 K too. Existence of 
the SSS below H c2 was predicted by H. Fink in 1965 [H- 
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FIG. 10: (Color online) Field dependence of x'l an d Xi f° r 
various temperatures. 
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FIG. 11: (Color online) Field dependence of X3 for different 
temperatures. 



face, the plane x — is in the center of the slab, and 
the external magnetic field is directed along the z-axis. 
Looking at the dimensionless order parameter in the form 
^(x, y, t) = 4>(x, t) exp(ifcy) the GL equation can be writ- 
ten as: 



HT c /T){-4> + |0M - + (a - kf^ = 0, Eq2 (2) 



da 
dx 2 



"2 \<f>\ (o-fe). 



(3) 



Here a is a y-component of the dimensionless vector po- 
tential. The order parameter is normalized with respect 
to the absolute value of the order parameter in zero field, 
the distances with respect to the coherence length at 
zero temperature, £oj ( x — * x /£,o, V — ► y/£o, I = 
L/£o) and the vector potential with respect to hc/2e^ 
(a = Aj (tic /2e£o))- The boundary conditions for calcu- 
lation of surface states are 0(0, t) = d(j>(±l,t)/dx = 
and a(0,t) = 0, da(±l,t)/dx = f)(i), where f)(t) is the 
dimensionless applied magnetic field. 

These nonlinear equations can be solved by numeri- 
cal methods. We add the time derivative d<j)/dt into 
the right side of Eq. ([2]) and seek the stationary solu- 
tions of the Eqs. ((2j[3|). Replacing the space derivatives 
by finite differences on the grid with step dx = l/N 
Eq. @ transforms into N first order differential equa- 
tions. The solution of the obtained linear algebraic sys- 
tem can be found by regular method. The grid with 
N=1000 points was used. In the surface state the order 
parameter differs from zero only near the surface, at a 
scale of several coherence lengths, £(T). Actually, the 
choice L = 5£(T) = D provides good accuracy for calcu- 
lating <j). The real dimensions of the investigated samples, 
L, considerably exceed this scale by 3-5 orders of magni- 
tude. Parameter k is not a gauge invariant quantity and 
we choose it using conditions a — at x — 0. In SSS the 
magnetic field in the bulk is constant. So we can obtain 
k for a thick slab with L ^> 5£(T) from the solution of 
the problem for a thin slab with D > 5£(T) by gauge 
transformation 



Increasing the DC field we can reach the field at which 
Xi or X3 becomes zero. This field can be considered as 
the third critical magnetic field H c ^. Both conditions ac- 
tually give the same value of H C 3. The experiment shows 
that the H c ^/ H C 2 ratio decreases with temperature. 



IV. THEORETICAL MODEL 

Let us consider a superconducting slab of thickness 2L 
in the parallel to its surface external DC and ac mag- 
netic fields. Due to the considerably short relaxation 
time of the order parameter [U [HJ one can use the sta- 
tionary GL equations. We choose the coordinate sys- 
tem in which the x-eods is perpendicular to the slab sur- 



k = k s + \) zs x (I - d) (4) 

and vector potential in the surface layer 

a(l — d + x) — a s (x) + \) zs x (I ~ d). (5) 

Here d = D/£q, index s corresponds to the problem for 
a thin slab, and \) zs is the z-component of magnetic field 
in the center of the thin slab. This note is important for 
numerical calculations. 



V. DISCUSSION 

It is well known (see, for example, [ill) that for a given 
external magnetic field there is a whole band of k for 
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which surface solutions exist. These solutions describe 
the noncquilibrium states and only one solution corre- 
sponds to the equilibrium state, for which the magnetic 
field inside the bulk equals its external value and the total 
surface current, J s , equals zero. Parameter k is an inte- 
gral constant of the nonstationary GL equations. That is, 
k is time independent, in contrast to <fi, in the frame of the 
GL model. The relaxation time of the order parameter <f> 
is considerably shorter than any ac period in our exper- 
iment. So when the external magnetic field is changing 
during the ac cycle, one may expect that <f> follows the in- 
stantaneous value of the magnetic field and k remains ap- 
proximately constant. Let assume that starting from an 
equilibrium state in some DC field, Hq, we increase the 
external magnetic field but simultaneously hold k con- 
stant. In this case the surface current J s becomes differ- 
ent from zero. It is possible to consider two definitions of 
the surface critical current J s i and J s i [l9, 34]. The first 
definition of such a critical current is J s \ = (c/4n)dh 3 i, 
where dh s i = H \ — Hq and Hi is the field for which the 
energy of the surface superconducting state equals the 
energy of the normal state [n| . The second definition is 
Js2 = {c/4ir)dh s 2, where dh S 2 = Hi — Hq and Hi is the 
field for which SSS disappears. The quantities dh s \ and 
dh s i have different values and different dependencies on 
the thickness of the sample, L. While dh s i dramatically 
depends on L, dh s i for L > 1000£ actually does not. The 
value of dh s i is considerably larger than dh s \ for large L. 
This difference is due to the large contribution of the 
magnetic field to the system energy, if the magnetic field 
in the bulk differs from the external field. These fea- 
tures are shown in Figs. 112a . 112b . where dh s i and dh s i 
are presented as a function of the DC magnetic field for 
different L's at T/T c = 0.9. In the reduced variables 
dh s i/H c i , dh s i/H c i, Hq/H c i the curves form is actually 
temperature independent. The assumption of slow relax- 
ing fc, permits one to understand qualitatively the effect 
of complete screening of a weak ac field with amplitude 
ho <C Hq in SSS. Ac surface current J s (k, H) is a function 
of the instantaneous values of the external magnetic field 
and k. This function can be calculated for a thin slab 
of several coherence length thickness and then using the 
gauge transformation, Eqs.([4]and[5]), to get a solution for 
a thick slab. As a function of k s and H, the J s (k s ,H) is a 
slow function of H. For example, at T = 0.9T C numerical 
calculation gives g^a^pS = 0.88 + 0.19(H ZS /H c2 - 1). 
Where H zs is magnetic field in the center of the slab. 
So for a thin slab, an almost complete penetration of 
the ac field inside the bulk takes place and the value 
of the surface current is very small. For a thick slab 
k =/= k s and the requirement of constant k during the 
ac cycle, implicitly means that k s also changed accord- 
ing to Eqs. ([4] and This leads to considerably large 
surface currents and to screening of the ac field. In real- 
ity, we have large dimensionless parameter L/£(T) that 
increases ac field screening. Fig. [13] demonstrates the cal- 
culated (in the assumption of constant k) \' = AM/dh s i, 
as a function of the DC field when the external field was 
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FIG. 12: (Color online) Field dependence of the surface criti- 
cal magnetic field (a) - dh a \{Ho/ H C 2) and (b) - dh S 2(Ho/ H C 2) 
for different slab thickness, L, at T/T c = 0.9 (see text). 



increased by dh s \. It is evident that for any macroscop- 
ically large sample, L > 5000£, the complete screening, 
X 1 = — 1/47T, should be obtained for DC fields excluding 
fields close to H C 3. However our experiments (Fig. 2]) 
do not confirm this conclusion. We see that x'i m f ne 
field H w (H c2 + H c3 )/2 already differs from -1/4tt. It 
means that slow relaxation of k takes place which leads 
to the losses and incomplete screening. 

For a given ac amplitude, x'i nas a maximum at some 
values of the DC field defined as H m (see Fig. [5]). H m was 
considered in Ref. [3l| as the DC field at which the ampli- 
tude of the ac surface current Jq s = (c/47t)ho equals ap- 
proximately to the critical value J s i. In order to test this 
in Fig. [T3] we show Jq s as a function of Hq and calculate 
a critical current J s i for a slab of thickness L = 5x 10 5 £. 
Theoretical data of the J s i were arbitrarily normalized in 
order obtain the intersection with the experimental curve 
at Hq/H c i = 1.25. While the theoretical dependence of 
J s i is almost a linear function of Hq, the experimental 
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FIG. 13: (Color online) Field dependence of x' f° r different 
slab thickness, L, at T/T c = 0.9 (see text). 



FIG. 15: (Color online) Out-of-phase susceptibility at max- 
imum, Xm, a s function of an excitation amplitude, ho- Inset 
shows the expanded view for weak ho in linear scale. 




H o' H c 2 



FIG. 14: (Color online) Surface critical current J s , found with 
assumption of Ref. [3l(. experimental values and calculated 
Jsi as a function of reduced DC field Hq/H C 2 (see text). 



curve starts from H /H c2 = 1.45 and is a nonlinear func- 
tion of H . One can conclude that losses observed in our 
experiment are not connected to the condition ho w dh s \ 
for Hq > H c2 . 

The maximal losses, x'mi a ^ H m , as a function of ho 
is shown at Fig. [TSl Inset to Fig. Qj}] shows that in the 
limit ho — > the losses do not disappear. In a linear 
system x" should be amplitude independent. While our 
experiments show a linear dependence x'm on the ac am- 
plitude (Fig. [T5]). It does not permit us to consider the 
response as a linear one even at very low amplitudes of 
excitation. Therefore more experimental measurements 
at low ac fields are needed. 

In general, the magnetic moment can be presented by 



following expression: 

M(t)= I K(t-t',h(t'))h(t')dt'. (6) 

J — oo 

For ho — 0.02 Oe the susceptibilities at higher harmonics 
are small and we can rewrite Eq. ([6]) as 

M(t)= [ K(t-t',h )h(t')dt' (7) 

J — oo 

considering only the response at the fundamental fre- 
quency. Under this approximation, the response at fun- 
damental frequency, matches the Kramers-Kronig rela- 
tions (KKR): 

x'i = + r f*' m 2 , dt (8) 

Jo ir{C, 2 - u 2 ) 

and then 

x"Mr°^#S)+xoo+ (9) 

poo 2u> 2n X "{0 J/- 

where u>o and u> m are the minimal and maximal avail- 
able frequencies in our experiment, respectively, and 
< zu < ujo- can be calculated from the available 

experimental data and be presented in the form 

"ma g 

I(u) = a + 61n|l -u 2 o/lu 2 \ + c n uj 2n . (10) 

n=l 

With c„ > we obtain x'i ( ro ) = ^h. Coefficients a, b and 
c„ could be found by least square fit. For ui 2 /ojf n << 1 
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FIG. 16: (Color online) Field dependencies of and pa- 

rameter a of Eq. dTOj) at T = 4.5 K and ho = 0.02 Oe (see 
text). 



it is sufficient to take into account only a few terms 
in Eq. (I10|) . Results of this approach are presented at 
Fig. [HI where Xi( w ) an d a as a function of the DC field 
are shown. The measured data in the frequency range 
15-1460 Hz xi at T = 4.5 K, was used for the calculation 
of I(oj) for 25 < w/2?r < 200 Hz with w /2tt = 17.5 Hz 
and w m /27r = 1455 Hz. The approximation of I(u>) by 
using expression Eq. (|10[) with n max = 1 produces Xii^) 
curve shown in Fig. 1161 Because c n tu 2 c n -i/uj m with 
uj 2 /uj'^ 1 w 0.02 one could expect that expression with 
ri raM = 1 gives the correct result. Taking into consider- 
ation the term C2 gives unphysical result, because C2 is 
very small and negative. It is due to the scattering of 
the experimental data and ignores in Eq. (|10p the depen- 
dence of zu on uj. Fig. [1751 shows that the calculated loss 
peak is approximately 3 times larger than the measured 
losses at oj/2tt > 20 Hz, Fig. 2] Qualitatively this be- 
havior can be explained as follows. Because x'{ exhibits 
a weak frequency dispersion we can estimate integral in 
the left side of Eq. © by 



R 



f 

- XiK + w /2) ln(^%5 



(11) 



In Fig. [17] we showed the correspondence between the 
estimated R by Eq. (fTTj) and the result of the numerical 
calculation of the integral in Eq. (jlip at Ho = 2000 Oe) . 
It is important that R has a positive sign for ui/2ir < 1000 
Hz and in the left side of Eq. (J9j) one has the sum of two 
negative values. So wc should expect a large contribution 
into the integral in Eq. |(5J) from frequencies outside the 
(luq, uj m ) region and the presentation of this contribution 
in the form of Eq. © gives a large value for the term b 
in Eq. dTOJ). 
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FIG. 17: (Color online) Frequency dependence of R, numeri- 
cal calculation and approximation by Eq. pip (see text). 



We believe that the observed in SSS losses are the re- 
sult of the relaxation k to its equilibrium value. This 
model can ascribe both the partial screening and losses 
for Ho > H C 2 . The other model assumes that the motion 
of the of 2D-vortices in the surface sheath [32j is responsi- 
ble for the losses [33[ . These vortices with surface density 
n s = Ho sin(0)/<^o appear if the applied field has a nor- 
mal component to the sample surface H n = f/osin(#), 
due to misalignment, or alternatively if the surface is 
not sufficiently smooth. One can estimate the conduc- 
tivity of the surface layer a = cr n H C 2/ Ho sm(8) where a n 
is the conductivity in the normal state. In our sample, 
a n w 10 17 CGS and the skin depth in the surface layer 
at frequency uj /2tt = 10 Hz is considerably larger for any 
real angle( ~ 10 -2 rad) to provide sufficient screening of 



the ac field by a layer with thickness 10 



io- 



The ac response of SSS resembles that of the spin-glass 
systems. Real and imaginary parts of xi can be well 
represented by a polynomial of In (a;) shown in Fig. llSl for 
Ho = 2 kOe and T=4.5 K. In this figure, presentation of 
X\i and xi ff by polynomial 



ao + ax ln(cj) + a2 ln 2 (w) 



(12) 



are shown for a considerably wide frequency region 15 < 
uj/2it < 1465 Hz. For some DC fields the coefficient 0,2 is 
small and one can get the spin-glass like xi ■ But x'{ also 
exhibits the frequency dispersion that is not typical for 
spin-glass systems. The "tt/2" rule H, x'l = -f Tf^jy, 
is not fulfilled in our data, Fig [19] 

The simple relaxation models of ac response is appli- 
cable only for a DC field near H c2 16]. If in analogy 
with a spin-glass system we assume that the magnetiza- 
tion moment of the sample, M(t), can be found from the 
relaxation equation: 



dM/dt = -vM - dh/dt, 



(13) 
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FIG. 18: Frequency dependence of x((y>) (panel (a)) and 
xiM (panel (b)) for H = 2 kOe (H = 1.25# c2 )at 4.5 
K. Continuous lines present the fit to the second order poly- 



nomial of ln( 


w) (Eq.(ll)). 




0.020 








Q 




0.015 






0.010 




-7t/2 dx'j/dln((B) 


0.005 


~& 






A 






A 





200 400 600 800 1000 1200 1400 
Frequency (Hz) 



FIG. 19: (Color online) The test of the "tt/2" rule for the 
experimental data at T = 4.5 K and Ho — 2 kOe. 



with subsequent averaging over the relaxation rates, then 



Xi = 



P{v)- 



-dv, 



(14) 



where P(y) is the distribution function of the relaxation 
rates. Using l/(u—iu) = J exp(— (i/— iw)i)dt we trans- 
form Eq.lO to 



I P{v) exp(— vt)dv, 
Jo 



(15) 



where P(t) = / °° P{v) exp(-vt)dv. So, if Eq.(3J 
describes adequately the experimental data with some 



FIG. 20: The test of the Eq. (TTo]) on the experimental data 
at T = 4.5 K and H = 2 kOe. 



P{v), then these two integrals should be equal each other 



P(t) = 2j ) °xi(w)cos(^)dw/7rw 
— 2 J Q Xii^sini^tydu/iruj. 



(16) 



Experimental data are available only for a finite fre- 
quency region 15 Hz < uj/2tt < 1465 Hz, while inte- 
grals in Eq. flT))) are expanded for all frequencies and we 
have to extrapolate our data to the entire frequency axis. 
This was done assuming that for lu/2tt < 15 Hz and 
uj/2tt > 1465 Hz Xi" is a power function of frequency ui p . 
As a result, the sin- and cos- Fourier transformations in 
Eq. p^)) give different values for P(t) as shown at Fig. [JD] 
where the ac response in H = 1.25H c2 was used. 

It is readily seen that the experimental data exclude 
the possibility consider the SSS as an analog of a spin- 
glass system. Equation ([7]) shows that in the quasilincar 
approximation the magnetization of the samples satisfied 
an integral equation 



G(t-t',h )M(t')dt' = h(t). 



(17) 



It is interesting to notice that the nuclear G(t, ho) can 
be extracted by the Fourier transformation of l/xi(w). 
Performing the same procedure as above, we obtained 
that sin- and cos- Fourier transformations in Eq. (fTTj) yield 
different values for G(t, ho) which is certainly due to 
the lack of experimental data for whole frequency axis. 
The extrapolation of the imaginary part of l/xi(w) gives 
more accurate results and we consider only the G(t, ho) 
that is obtained by the sin-Fourier transformation of 
1/Xi(w). Good approximation of G(t,ho) provides the 
expression G(t,ho) — A(t)/t q with slow function A(t) 
for t > 7r/1465 = t c . For t < t c function G(t,h ) is 
singular, but integral J Q C G(t, ho)dt has a finite value. 
The parameters q and A(t), depend on the DC field. 
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FIG. 21: Time dependence of nuclear of Eq. (|17[) for several 
DC fields near the absorption maximum at T = 4.5 K. 



For example, in field Hq = 1.25H C 2 q = 0.876 and 
A(t) = - exp(l. 285 - 0.00842 ln 2 (27rf)). In Fig. [21] we 
show G(t, ho) for some values of the DC magnetic field 
and the inset presents q versus H . So, the dynamics 
of SSS is governed by an integral equation with retarda- 
tion. This feature distinguishes SSS from other known 
systems. 

VI. CONCLUSION 

In this paper we have studied the low frequency lin- 
ear and nonlinear dynamics of the SSS of a single crys- 
tal of yttrium hexaboride. The tunneling spectra were 



studied as well. Tunnel measurements allow us to make 
the assumption, that in this single crystal, unlike ZrBi2, 
near the surface the electron-phonon interaction is sup- 
pressed and the situation of weak coupling is realized. 
We showed that the surface superconducting states de- 
fine the peculiarities of the low frequency response. In 
spite of different behavior under magnetic fields (ZrBi2 
is a type-I superconductor and YBq is that of type-II) 
and different surface properties the two materials exhibit 
very similar and universal ac characteristics reflecting the 
nature of the SSS. In both cases we observed a nonlin- 
ear response for very weak ac amplitudes (in experiments 
with YB 6 ho was as small as 0.005 Oe) and the ques- 
tion about the existence of a linear response is open. An 
extrapolation of the low-amplitude data did not reveal 
a linear regime. Similar to spin-glass systems (where 
finite losses at considerably low frequencies exist), the 
real part of the susceptibility exhibits a logarithmic fre- 
quency dependence at some DC magnetic field. But the 
out-of-phase component has a frequency dispersion. The 
frequency dispersion in SSS is different from that of the 
spin-glass systems. The slow relaxation of the phase of an 
order parameter leads to a frequency dispersion of the ac 
susceptibility. The analysis of the experimental data by 
means of Kramers-Kronig relations allow us to make the 
assumption of the presence of the loss peak at frequencies 
below 5 Hz. 
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